Ergodic transitions in continuous-time random walks 
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We consider continuous-time random walk models described by arbitrary sojourn time probability 
density functions. We find a general expression for the distribution of time-averaged observables for 
such systems, generalizing some recent results presented in the literature. For the case where sojourn 
times are identically distributed independent random variables, our results shed some light on the 
recently proposed transitions between ergodic and weakly nonergodic regimes. On the other hand, 
for the case of non-identical trapping time densities over the lattice points, the distribution of time- 
averaged observables reveals that such systems are typically nonergodic, in agreement with some 
recent experimental evidences on the statistics of blinking quantum dots. Some explicit examples 
are considered in detail. Our results are independent of the lattice topology and dimensionality. 
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I. INTRODUCTION 

Boltzmann's hypothesis of ergodicity is a central con- 
cept in Statistical Mechanics. Roughly speaking, in an 
ergodic system, for a long time observation, the resi- 
dence time of a trajectory in a given region of the phase 
space is proportional to the volume measure of the re- 
gion. Despite the high success of Boltzmann's descrip- 
tion of large systems, the ergodic hypothesis cannot be 
used, for instance, for systems whose phase space can 
be subdivided in mutually inaccessible regions. A sub- 
tler physical scenario of nonergodicity was introduced by 
Bouchaud pj] in the context of glass dynamics: the so- 
called weakly nonergodic systems are also nonergodic, but 
their phase spaces are not subdivided in mutually inac- 
cessible regions. Recently, there has been great interest 
in particular, in the weak ergodicity breaking phe- 
nomenon, where a transition to an ergodic phase may 
occur. Notice that Bouchaud's ideas about weak ergod- 
icity breaking are, in turn, closely related to concepts 
that have been previously considered in the mathemat- 
ical literature about ergodic theory and stochastic pro- 
cesses, see, for instance Q. In particular, a weakly non- 
ergodic regime corresponds to a situation in which the 
state space of a (semi)-Markov process is connected, i.e., 
any state can be reached from any other state with finite 
probability in a finite number of steps, but the fraction 
of occupation time for a given state is not equal to its 
invariant spatial measure. 

From the physical point of view, weakly nonergodic 
systems have proved to be relevant in many applications 
as, for instance, complex networks j|, weak turbulence 
01 , and, in particular, they are at the basis for the statis- 
tical modeling of atoms trapping by laser cooling devices 
Q. In these models, the atom quantum dynamics 



are equivalent to a classical random walk in the momen- 
tum space, where the standard deviation of the jump 
lenghts Ap is of the same order of the incident photon 
momentum. The trapping process consists effectively in 
successive frontal collisions, for some given time interval, 
between the atom and the resonant laser photons, in a 
process called subrecoil cooling. Statistically, it is equiv- 
alent to the application of a controllable potential jump 
rate R{p), responsible for the trapping near p = 0. The 
resulting scenario is a kind of continuous-time random 
walk where rare events play a dominant role, with Levy 
type non-gaussian probability density functions (PDF) 
tp{T) ^ 74Qr^(^+"-' governing the trapping times. The 
divergence of the average trapping time 
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for < a < 1, with E ( ) denoting the expectation value 
with respect to the PDF ipir), is precisely at the root 
of the subrecoil cooling mechanism effectiveness. Such 
mechanisms are typically nonergodic, and the origin of 
the nonergodicity is usually attributed to the divergence 
of the average trapping time \§\ . 

We investigate here the PDF of time-averaged observ- 
ables for a general class of continuous-time random walk 
models (CTRW), extending considerably the classes al- 
ready considered in [2]. The nonergodic properties of 
CTRW models are studied here by comparing PDF of 
ensemble averages with fiuctuations of time averages. Ini- 
tially, we consider CTRW models described by an arbi- 
trary trapping time PDF 'iJj{t), in which all the sojourn 
times r of the L lattice points 1 < x < L are identically 
distributed independent random variables. We show that 
models with finite f are always ergodic. In the cases 
where f diverges mildly, one still has ergodicity, but for 
stronger divergences the dynamical regime is weakly non- 
ergodic and characterized by a Lamperti type PDF of 
time-averaged observables, in agreement with the recent 
results announced in [2| . We also extend our approach to 
include CTRW models with non identical trapping time 
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PDF over the lattice points. For such cases, we show that 
even models with f finite for all lattice points are typ- 
ically nonergodic. The obtained PDF of time-averaged 
observables clarify some nonergodic features reported re- 
cently for blinking quantum dots systems [lol . |Tlj . 

The paper is organized as follows. In Section we 
provide the basic concepts and definitions required to 
describe ergodic and nonergodic properties of CTRW 
models. In Section IIIIl the distribution density of time- 
averaged observables for the CTRW model with identical 
trapping time PDF is obtained. Section IIVI is devoted 
to discuss the weak ergodicity breaking phenomenon, ex- 
tending some of the results obtained in [3| . Some explicit 
examples of PDF are considered in Section [V] In Section 
rvTl we investigate the CTRW with non identical trap- 
ping time PDF. We conclude in Section I VIII with some 
remarks about possible universal aspect of nonergodic 
fiuctuations close to the weak ergodicity breaking. 



where 



[0] = 



(4) 



with the square bracket denoting, hereafter, the average 
of a given operator with respect to the set . . . 
The generating function ([3]) is related to the correspond- 
ing density function of time average observable O as 



p{d) = E{s{6-[0])) 
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Notice that <^ must always be chosen in order to assure 
that l3 + ^[O] > for any set {ii, . . . , ii}. 



II. STATISTICAL APPROACH FOR CTRW 
MODELS 

Here, we consider a general CTRW model on a lat- 
tice with points x = 1, . . . , L, subject to certain trapping 
statistics. The topology and the dimensionality of the 
lattice are irrelevant for our purposes. The notion of 
random waiting time between successive steps was origi- 
nally introduced by Montrol and Weiss [l^] ■ In the model 
discussed here, the particle can jump to one of its nearest 
neighbors after waiting for a random time r. The trap- 
ping statistics is given by the waiting time —J2^=i'''i^ 
where Ti is the «th sojourn time of a given lattice site x. 
The set {ri, . . . , r„^} for 1 < x < L is composed of non- 
negative and identically distributed independent random 
variables with a common and arbitrary PDF iP{t). As 
in Ref. we also consider here = nP^'^ for the 
equilibrium regime of the random walk, which is reached 
after a sufficiently large number of jumps n. The Laplace 
transform for the random variable tx is given by 
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%l){Ti)e 
exp {-nP^^^hiu)) . 



(2) 



Evidently, h{u) is a monotonically increasing function for 
M > 0, with h{Q) = 0. 

In order to investigate the ergodic phase transitions 
for the trapping process, it is interesting to introduce 
the generating function p(^) developed in Ref. and 
successfully used in Ref. ^ for describing time aver- 
age observables. We will assume here, however, that the 
generating function depends on an additional parameter 
/3 > as follows 



p(/3,0=i^((/3 + ^[O])- 



(3) 



A. Ergodicity and weak nonergodicity 

Physically, we introduce the idea of a weakly noner- 
godic behavior by demanding that the lattice be not sub- 
divided, from the dynamical point of view, in mutually 
inaccessible regions. This is equivalent to impose that [3] 



< Pf * < 1 for all 1< X < L. 



(6) 



This condition is common to ergodic and weakly noner- 
godic systems. In fact, we say that a CTRW is ergodic if 
the condition ([6]) holds and, besides, the PDF ([5]) could 
be written as p{6) = pe{0) = 5{6 - {O}), where 



(7) 



will denote, hereafter, the equilibrium ensemble average, 
being Oj the value of the observable O when the particle 
is at the lattice point j. The PDF for the ergodic regime 
can be rewritten as 



PEiO) 
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Combining Eqs. ([5|) and ^ we obtain the corresponding 
ergodic generating function: 



1 



(9) 



The generating function ([9|) will be used hereafter as an 
unequivocal mark of ergodicity. In fact, we call weakly 
nonergodic a system which obeys (jS]) but for which the 
generating function ([5]) is not equivalent to (|9]). In other 
words, a weakly nonergodic system is a system where it 
is possible to get to any state from any state, but still 
time and ensemble averages do not coincide. 
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III. DISTRIBUTION DENSITIES OF TIME 
AVERAGED OBSERVABLES 

The generating function defined by Eq ^ can be 
rewritten in the following form 



where 



/5(/3,C) 



pOO pOO pOO pOO 

/ ds dt dti-^itx)... I dtL'i'itL) 
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5 t-^t, 



(10) 



where ^{tx) is the corresponding PDF for t^- We also 
demand 



Uj= (3 + iOj > 0, l<j<L, 



(11) 



in order to guarantee that /? + ^[O] > 0. Now, by using 
the Laplace transformation given by Eq. ([2]) and taking 
into account that 




1 

2^ 



dk e 



(12) 



Eq. (ITUl) can be cast, after successive changes of variables 
similar to those one performed in Ref. jH, in the form 



OO /"OO 



dsdk ■ 



2Ti{k-if3sY 



(13) 



The fc-integration can be performed by using the Cauchy 
formula, leading finally to 



m = - r d.s-4-e'^^^'^<^^''^l 
Jo s dp 



(14) 



Equation (ITi|) is similar to that one found in Ref. [2j 
for the particular, but very important, case of one-sided 
Levy PDF, for which h{u) — Cu", with < a < 1. Here, 
however, we did not use either the generalized central 
limit theorem or scaling hypothesis. 

In order to solve equation (|14p for arbitrary functions 
h(u), we introduce a new parametrization for the integral 

fi{s)^{h{sU)), (15) 

from which it is always possible to obtain s = s(/i) since 



{Uh'{sU)) > 0. 



(16) 



Notice that /i(0) = and /i(oo) = fimax- With the new 
parametrization. Equation (|14p reads 



(17) 
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/(m) = 



{Uh'{s{n)U))' 



(18) 



The integral pT)) can be evaluated in the limit n oo, 
leading to (see Appendix \K\ for details) 

m = hm /(/.) = lim (19) 

from which the PDF of time averaged observables p{0) 
can be obtained by using Eq. ([5]), generalizing the previ- 
ous results obtained in [2| . Notice that Eq. depends 
only on the behavior of h(u) near u = 0+, confirming 
that only the asymptotic tail behavior of iP{t) is rele- 
vant in the limit n oo. It should also be emphasized 
that the sojourn time r is a non-negative unbounded ran- 
dom variable. Hence, only functions exp(— /i(u)) that are 
Laplace transforms of PDF V(t) with support on [0 oo) 
are relevant here. According to Bernstein's theorem [14| . 
the functions exp(~h{u)) must be completely monotonic, 
i.e., they should obey 



(-l)'=^exp(-M^^))>0, 



M > 0. 



(20) 



IV. 



ERGODIC AND WEAKLY NONERGODIC 
REGIMES 



The first conclusion that one can draw from Eq. 
is that systems with finite average trapping time f are 
ergodic since for such cases lim„_j.Q+ h'(u) = f, imply- 
ing that Eq. (fT9|) reduces to Eq. ([9|). On the other 
hand, the relevant weakly nonergodic models with di- 
verging average trapping time are those ones considered 
in [2|, for which h{u) = Cu", for small non negative u, 
with < a < 1. The associate PDF in this case are the 
well known stable Levy densities, for which the asymp- 
totic behavior for large r is given by 



V'(T)^A„r-(i+"). 
For such models, Eq. reads simply 



m = 



(21) 



(22) 



leading to a density function of Lamperti type [2]] , where 
the ergodic regime can be recovered in the limit a — >■ 
1~. Interestingly enough, many other subtler ergodic 
solutions do also exist. This is the case, for instance, 
of the function h(u) = — Culogu for small u, for which 
Eq. p9|) also reduces to Eq. ([9]). This case corresponds 
namely to a PDF that asymptotically tends to the Levy 
density pT|) with a — \. 

Since PDF with finite average sojourn time f give rise 
to ergodic behavior, it would be worthy to classify the 
possible PDF t\}{T) with diverging f in order to identify 
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possible nonergodic regimes. Since one must demand 
E{1) ~ 1, all PDF ?A(t) shall asymptotically decrease 
faster than r^^. On the other hand, in order to have 
f — E{t) diverging, il^ij) cannot decrease faster than 
T^^. Hence, a PDF with diverging average sojom'n 
time f must obey 



— < V(t) < — , 
r r 



(23) 



for large r, with Aq and Ai arbitrary positive constants. 
The cental point here is that the class of power-like func- 
tions like h{u) = Cw", with < a < 1, which correspond 
to the Levy distributions (f2T|) . does not exhaust all the 
interval (j23p . In particular, the lower bound of the inter- 
val does not correspond to any of these functions. As we 
have already shown, such case (ergodic and with a = 1) 
corresponds indeed to the function h(u) = —Culogu for 
small u. PDF with asymptotic behavior of the type 



V'(r) 



r log'' 



(24) 



where i?^ and 7 > 1 are constants, obey (f23| . They, in 
fact, accumulate in the upper bound. On the other hand, 
PDF of the type 



■0(t) 



(25) 



where d, and 1/ > are constants, also belong to the 
interval and accumulate in the lower bound. 

In order to extend our analysis to consider the PDF 
with logarithmic terms are those ones of ([M)) and ([25]) . 
one can make use of Karamata's Abelian and Tauberian 
theorems [31 for the Laplace-Stieltjes transform 



e-"^dT (r) , 



for M > 0, which states that 



-h(u) 



r(p + i)T - 

' u 



for u — ^ 0-I-, if 



T->-00 1 [t) 



(26) 



(27) 



(28) 



In the present case, since T (r) is the cumulative distri- 
bution function associated to tpir), the condition ([28)) is 
automatically fulfilled with /? = 0. Hence, from (P?]). we 
have 



h{u) 



1-T ( - 

u 



(29) 



for {)+. 

For the PDF dM]), we have, by using ^ 



h{u) 



1 



7 - 1 I logu|T-l ' 



(30) 



for small u and 7 > 1. For this case, we have h'{u) 
B^/{u\ logul'') for small u, implying, from ()19|) . that 



m = {u- 



(31) 



for any value of 7 > 1, which, interestingly, coincides 
with the limit a — > of the Lampertian case given by 
(|22l) . For the PDF ([25l), on the other hand, we have 



h{u) ^ Ci,u\ loguj''. 



(32) 



for small u and j/ > 0, which give rises to a ergodic 
generating function 



p{u) = {uy 



(33) 



for any > 0, which, of course, also coincides with the 
limit a — >■ 1 of the Lampertian case (|22p . The Lamper- 
tian generating function ([2^ . with < a < 1, seems to 
be enough to describe any CTRW of the type considered 
up to here. 



V. STABLE PDF 

Any physical application of the preceding sections re- 
sults would require, of course, stable PDF '0(r). We 
recall that a stable PDF S'q(t; cr, /3, /i) is characterized 
by four parameters: < a < 2, which determines the 
asymptotic falling tails; tr > is the corresponding scale; 
— 1 < /3 < 1 being the skewness; and ji the shift param- 
eter. The support of a generic stable density is given by 

M 



[n,oo), < a < = 1, 
(—00,^], < a < l,/3 = 
(— oo,+C)o), otherwise. 



-1, 



(34) 



Since r is an unbounded non-negative random variable 
we disregard the case /3 = — 1. The natural stable PDF 
in this case are 



V^t) = 5a(T;cr, 1,0), 



(35) 



with < a < 1, which have the asymptotic behavior 
given by (|21|) . There are, however, several other useful 
possibilities. 

Let us redefine the support of a stable PDF by taking, 
first, /i = C > and then t — > |t— C| so that, for (3 ^ — 1, 
one has 



^(r) 



iV5'a(|r-C|;cr,/3,0), for t > 0, 
0, for T < 0, 



(36) 



where N is the pertinent normalization constant and < 
a <2. The general asymptotic behavior of the PDF (|36|) 
is given by (I^Tt 1?^]. The average sojourn time f — E{t) 
diverges as t^~" for < a < 1, and logarithmically for 
a = 1. According to the results of the previous section, 
the cases for which f is finite are ergodic. Some explicit 
examples will help to illustrate our main results. 
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A. A non-negative delta sequence PDF 



where 



Consider the non-negative delta sequence PDF on 
[0, cx)) given by 



V'(t) = N 



sin^ 7(t - C) 
7(r-C)2 ■ 



(37) 



where TV is the appropriate normahzation constant and 
7 and C are positive parameters. PDF ([57| is the one- 
sided version of the density defined in Ref. In the 
hmit 7 — >■ oo, it corresponds to a delta function centered 
in r = C. The function h{u) associated to (|57)) is given, 
for small non-negative u, by 



D = 2N 




^ -f (1 - N)a 



(43) 



For large times, the term (cru)^/^ prevails, leading to the 
{a = 1/2) Lamperti's statistics (f22|) . This result coin- 
cides with that one obtained in Ref. 2] for the PDF 
([55)1. Note, however, that cr — > leads to an ergodic 
regime, regardless of having a — 1/2! This is a some- 
what surprising result. 



VI. NON-IDENTICAL TRAPPING TIME PDF 



where D is a constant depending on 7 and C. Despite 
of having a diverging first moment f , the PDF ([37]) gives 
rise to an ergodic regime since the dominant term in the 
function h{u) for small u is precisely —ulogu, that it is 
known from the results of the last section to reduce ^T9\\ 
to the ergodic generating function (|9]). 



B. The Cauchy stable PDF 



The Cauchy stable PDF, defined for r e [0, 00) as 



V'(r) = N 



(r - C)2 + a2 



(39) 



where N is the appropriate normalization constant and 
a and C are positive parameters, is another example of 
a PDF with diverging f but that, nevertheless, gives rise 
to an ergodic regime. The function h{u) in this case is 
given by 

h{u) = {D - N\og{cru)) cru + 0{{auf), (40) 

for small and non-negative u, where 13 is a constant de- 
pending on C and a. Again, we have a situation where 
the logarithmic term dominates and renders (|19p to the 
ergodic expression 



C. The Levy stable PDF 

Based on the one-sided PDF let us consider the 
(a = 1/2) Levy density with support on [0, 00) given by 



■ exp 



■2\t-C\ 



C\3/2 



(41) 



where N is the appropriate normalization constant and a 
and C are positive parameters. The corresponding h(u) 
function can be calculated for u — > 0"'' as 



So far we have considered only situations where a sin- 
gle trapping time PDF tpiT) is sufiicient to describe the 
random dynamics. However, there are situations where 
the trapping times are not identically distributed over the 
lattice points. This is the case, for instance, of the fluo- 
rescence blinking observed in some colloidal nano crystals, 
namely the case of certain quantum dots systems jlo[.[l]|. 
Once a laser pulse incides on a quantum dot of these 
systems, its fluorescence intensity randomly switches be- 
tween bright (on) and dark (off) states. The blinking 
quantum dots are typically characterized by means of 
the statistics of on/off times, whose distributions exhibit 
power law decay ipk{T) ^ AkT~'^^~^°'''\ where k corre- 
sponds to the on/off states and < < 1. Nonergodic- 
ity for such systems has been reported from experimental 
observations [To| . 

Our approach can be generalized to include such kind 
of situation. To this end, let us consider now a CTRW 
model described by a set of arbitrary trapping time PDF 
iJjj{T) governing the sojourn times Tj of the L lattice 
points 1 < j < L. In such non-homogeneous lattice, the 
sojourn times are not identically distributed anymore, 
but they are still independent random variables. We can 
extend all the results Sections |ll] and IIIII essentially by 
replacing h{u) by hj{u), resulting finally in 



m 



lim 



yL pe 



0- Y.- Pr^,h){sU,) 



(44) 



h{u) = nV2^ +{C- D)u + 0((o-u)3/2), 



(42) 



instead of ([T9|) . The first conclusion we get from (|44|) is 
that weak nonergodicity is also present in this case, since 
we might have nonergodic regimes for which < Pj"^ < 1 
for all 1 < j < L, assuring that the lattice is not dynam- 
ically subdivided in mutually inaccessible regions. Fur- 
thermore, we have also that the weak ergodicity break- 
ing is not a structurally stable phenomenon since ergodic 
transitions are impossible for systems where the sojourn 
times T are not identically distributed. Even for systems 
where the average values fj are finite for all lattice points, 
but different, one has the predominance of nonergodicity. 
In order to illustrate this scenario, some explicit exam- 
ples are useful. For any of the three cases depicted in 
Table [H we have that the expression for p{0) as defined 
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hjju) ^j{t) Tj 

CjU faster than finite 
—Cjulogu AjT~'^ logr 



TABLE f : Some possible non-liomogeneous PDF over the lat- 
tice. The first column depicts the behavior of hj{u) for small 
non- negative u, corresponding to the lattice point j. The 
second column has the associate PDF asymptotic behavior 
for large t. The last column shows how the average sojourn 
time diverges. We assume that Cj are positive constants and 
< a < 1. 



ijki is infinite or vanishes, then, respectively, state k pre- 
vails over I, or vice versa. In the case where prevailing 
states exist, only them are relevant for the calculation of 
the density of time averages. Notice that even if a state k 
has a visitation fraction negligible when compared to an- 
other state I, i.e. P^' <C Pf'', state k can still prevail over 
state I if rjki = oo. Irrelevant states can be visited many 
times, but the time spent among them by the CTRW 
dynamics is negligible and, hence, they should not con- 
tribute to the time average. This is the mechanism be- 
hind the nonergodicity of non-homogeneous lattices. 



by (O are obtained from the corresponding single trap- 
ping time PDF cases by replacing PI 



,eq 



replacing 



(0)-> 



(CO) 
(C) 



Pf'^C,, I.e., by 



(45) 



what prevents any transition between weakly nonergodic 
and ergodic regimes. We can also conclude that the aris- 
ing of narrow PDF does not necessarily mean that a sys- 
tem has reached the ergodic regime. 

Another interesting dynamical aspect of the large time 
evolution of systems with non-homogeneous PDF is the 
competition between the lattice points. Let us introduce 
the survival rate rjki between two lattice points k and I, 
defined as 



rjki 



.... 

ti-i.o+ h'i{u) 



lim 



(46) 



The survival rate is, of course, related to the ratio ff^/fi 
for the lattices points k and I. Such relation, however, is 
not obvious for general PDF. Notice that 

/>oo 

e-''^^"'>h'^{u)^ / e-"^dA,(T), (47) 
Jo 



where 



Ai(r) = / sipi{s) ds. 



(48) 



Provided that Ai(r) obey a condition like (|28p. we have 
by using Karamata's theorem for (|T7)) 



1)A. - 
' u 



(49) 



for u — > 0+. Since Ai(oo) = f,;, we conclude that the 
survival rate (j46p does no coincide, in general, with the 
ratio fk/fi. Nevertheless, the survival rate is indeed 
an apropriate quantity to compare the residence time 
of the CRWT dynamics in different points of the non- 
homogeneous lattice. 

According to (j46l) . two types of behavior can occur. 
Both states coexist if r/hi is finite and non-negligible. If 



VII. CONCLUDING REMARKS 

Weak ergodicity breaking has been investigated 0- 
01 by considering Levy PDF (|2ip . for which the time- 
average densities are given by a Lamperti distribution by 
means of ([2^ . For this class of PDF, the weakly noner- 
godic regime breaks up, giving origin to an ergodic one, 
in the limit a — >■ 1. Despite that our results strongly in- 
dicate that Lamperti distributions seems to be enough to 
describe any time-average density for CTRW with PDF 
obeying (|23p . our results also show that the limit a — > 1 
is not the only way to obtain ergodic regimes. This fact is 
explicitly illustrated by the PDF given by (|4ip , for which 
an ergodic regime arises for cr — >■ 0, irrespective of having 
Q = 1/2 for such PDF. This results challenges the naive 
association between ergodicity/ weak nonergodicity and 
diffusion/subdiffusion as suggested in jj], which indeed 
appears reasonable at first sight since models with PDF 
(1211) typically exhibit anomalous diffusion characterized 
by (a;^) ~ t" . The characterization of any possible uni- 
versal behavior close to the weak ergodic breaking cannot 
by achieved by analyzing only the breaking associated 
with a — ^ 1. 

As to the case of non-identical sojourn times over the 
lattice, the issue of weak ergodicity breaking is still more 
involved. Even for PDF with finite average sojourn times 
fj, we do not have an ergodic regime, in contrast with 
the case of homogeneous lattices, for which a finite f will 
necessarily imply in an ergodic regime. This nonergodic 
regime can be considered also as weakly nonergodic since, 
from pS)) . we have that the phase space is similar to the 
homogeneous case being, in particular, not subdivided in 
mutually inaccessible regions. Ergodicity, in this case, is 
recovered in the limit of equal PDF. This can be illus- 
trated by a simple explicit example. Let us consider a 
Markov chain with only two states A and B and with 
the transition matrix 



(50) 



The corresponding invariant weights are — Pg'' = 
1/2. Assume now that the sojourn times and Tg are 
governed by the PDFs V'a(t) = A^e"^-*"^ and ^b{t) = 




-\bt 



After n ^ 1 steps, we will have ua ~ tib 
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n/2, but the visitation fraction of the dynamics in the 
two states will be in general different since 



nBE{TB) Xa' 



(51) 



For the case of diverging fj , besides of the consequences 
associated with (|45|) . we have also new features associated 
with the possible divergence or vanishing of the survival 
rate ([^5]) . This leads to the possibility of a new weakly 
nonergodic regime, for which, despite of the CTRW 
spreading over all lattice points, the time-averaged dis- 
tributions will depend only on some points, namely the 
prevailing ones according to the survival rate De- 
spite these points are not surprising from the probabilis- 
tic point of view, they certainly deserve a deeper physical 
investigation. 
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Appendix A 

The first observation about the limit of large n of equa- 
tion P?)) is that one cannot apply Watson's lemma [13] 
directly since we cannot assure, a priori, if the deriva- 
tives of /(^) are finite or not for /i — > 0+. Let us take, 
for instance, the first derivative 



hm /'(/.) 



( {Uh"{sU)) ,,,, ...AU^h^^imi 

(Al) 

where (fTH]) was used. Although, apparently, for all phys- 
ically relevant functions h{u) the limit (jAl|) is indeed 
finite, one cannot rule out, in principle, possible situa- 
tions where it might diverge. Fortunately, thanks to the 
boundedness of /(/i), one can evaluate (|17p without using 
Watson's lemma. 

By introducing the new variable v — nii, Equation (jl7p 



can be cast in the limit n — oo as 



p = lim 







' dv, 



(A2) 



with f{fi)>0 given by (ITSl) . First of all, let us suppose 
that ^max < oo- Equation (jA2p can be decomposed in 
this case as 



p= lim Io{n) + lim Ii{n), 



where 



and 



dv 



h{n) = I f(l)e--dv. 

For both integrals, we have 

fr9i{n) < I,{n) < ftg^in), 
with i = 0, 1, where 



9o{n) = 1 - e 
9i{n) 



where In = 



h = inf /(m), 
ft = sup /(m), 



0, 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
(A8) 
(A9) 

(AlO) 



Since 



and Ii = 

/(/i) is bounded, we have < < oo. Taking 
into account that lim„_>oo (n) = 0, one has from (|A6p 
that Imin^oa Ii(n) = 0. As to the integral Io{n)^ no- 
tice that lim„_>oo /cT = lini„^oo fo = Hm^-^o+ /(m) and 



lim„_>.oo 5o('^) = Ij implying finally from (|A6|) that 



p — lim n 



f(p)e-^^dp= lim, f(p). (All) 







The cases for which ^max diverges can be treated in an 
analogous way, by choosing \pnp^ with finite /x* , instead 
of y/nfj,rnax, for the decomposition of (|A2|) into Io{n) and 
Ii{n) integrals, leading to the same final resuh (|jlT|) . 
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